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For any realistic halo profile, the Galactic Center is predicted to be the brightest source of gamma- 
rays from dark matter annihilations. Due in large part to uncertainties associated with the dark 
matter distribution and astrophysical backgrounds, however, the most commonly applied constraints 
on the dark matter annihilation cross section have been derived from other regions, such as dwarf 
spheroidal galaxies. In this article, we study Fermi Gamma-Ray Space Telescope data from the 
direction of the inner Galaxy and derive robust yet stringent upper limits on the dark matter's 
annihilation cross section. Even for the very conservative case of a dark matter distribution with 
a significant (~kpc) constant-density core, normalized to the minimum density needed to accom- 
modate rotation curve and microlensing measurements, we find that the Galactic Center constraint 
is approximately as stringent as those derived from dwarf galaxies (which were derived under the 
assumption of an NFW distribution). For NFW or Einasto profiles (again, normalized to the mini- 
mum allowed density), the Galactic Center constraints are significantly stronger. For example, for a 
dark matter particle with an annihilation cross section of av — 3 x 10 -26 cm 3 /s to bb and distributed 
according to an Einasto profile, we robustly rule out dark matter particles with a mass of less than 
285 GeV. 

PACS numbers: 95.85.Pw,95.55.Ka,95.35.+d; FERMILAB-PUB-12-505-A 



I. INTRODUCTION 

A major scientific mission of the Fermi Gamma-Ray 
Space Telescope (FGST) is the search for gamma-rays 
from dark matter annihilations. To date, Fermi data has 
been used to search for dark matter annihilation products 
from dwarf spheroidal galaxies [U [2] , galaxy clusters |3l- 
[5], the Galactic Halo [B], galactic subhalos [7J, and from 
among the isotropic gamma-ray background 8, 9J. The 
most stringent of these constraints are beginning to probe 
annihilation cross sections at or around the value pre- 
dicted for a simple thermal relic (av ~ 3 x 10 -26 cm 3 /s) 1 , 
at least for dark matter particles with masses below a few 
tens of GeV and which annihilate to final states which 
result in significant fluxes of gamma-rays. And while the 
particle (or particles) that make up the dark matter may 
possess an annihilation cross section that is well below 
this value, the range of cross sections that is presently 
being probed by Fermi represents a very well motivated 
and theoretically significant benchmark. 

Observations of gamma-rays from the inner region of 
the Milky Way are of particular interest for indirect dark 
matter searches. Due to its relative proximity and the 



The actual value of the annihilation cross section required in or- 
der for a simple thermal relic to be produced with the measured 
dark matter abundance is actually slightly smaller than the com- 
monly used canonical value; approximately 2.2 X 10 -26 cm 3 /s for 
dark matter particles heavier than about 10 GeV |10| . 



high densities of dark matter present in this region, the 
Galactic Center is predicted to be the single brightest 
source of dark matter annihilation products in the sky. 
For dark matter distributed according to the commonly 
used NFW (Navarro-Frenk- White) profile, for example, 
one predicts a flux of gamma-ray annihilation products 
from the Galactic Center that is roughly four orders of 
magnitude higher than from even the most promising 
of the known dwarf spheroidal galaxies (Segue 1, Ursa 
Major, etc.). Even if the dark matter distribution in 
the inner Galaxy has a significant constant-density core, 
the Galactic Center will significantly outshine all other 
sources of dark matter annihilation products. Temper- 
ing this advantage is the fact that astrophysical back- 
grounds from the Inner Galaxy are fairly bright, and are 
not particularly well understood. So while even a very 
faint gamma-ray signal from one or more dwarf galaxies 
could possibly provide us with evidence of annihilating 
dark matter, a much higher flux of annihilation products 
would be required from the inner Galaxy before it could 
be identified as such. 

In previous studies, a number of independent groups 
have studied Fermi data from the direction of the Galac- 
tic Center |11H16| (preliminary results from the Fermi 
Collaboration have also been presented [T7]-(Tnj). In 
particular, in Refs. [ITJ [12j [16], a bright gamma-ray 
source was found to be present at the Galactic Center, 
with spectral and morphological characteristics consis- 
tent with those predicted from annihilating dark mat- 
ter. This conclusion was recently confirmed by Ref. [T3] . 
A debate regarding the origin of this emission is ongo- 
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ing and has focused on a number of possibilities, includ- 
ing that these gamma-rays largely result from cosmic ray 
collisions with gas [TP - THl |2"01 [2Tj , from a population of 
~ 10 3 millisecond pulsars [TTHT31 |2"2"1 12"3"] . or from anni- 
hilations of ^7-40 GeV dark matter particles with an 
annihilation cross section on the order of av ~ 10~ 26 
cm 3 /s [TiTfTBI [TB] . Throughout this paper, we will re- 
main agnostic as to the origin of this emission, and use 
the observed spectral and spatial distribution of the ob- 
served gamma-rays to derive upper limits on the dark 
matter annihilation cross section. We find that even for 
very conservative choices for the dark matter distribu- 
tion (such as cored halo profiles), wc derive constraints 
on the dark matter's annihilation cross section which are 
comparable to the strongest constraints found from other 
regions of the sky. In particular, our most conservative 
constraints from the Galactic Center (assuming a distri- 
bution with a kiloparsec-scale constant-density core, nor- 
malized to the minimum value compatible with the Milky 
Way's rotation curve) are comparable to the constraints 
derived by the Fermi collaboration from the combination 
of all dwarf spheroidal galaxies (under the assumption 
of NFW halo profiles). If instead we adopt an NFW, 
Einasto, or contracted profile to describe the dark mat- 
ter distribution in the Galactic Center, the constraints we 
derive on the dark matter annihilation cross section are 
significantly more stringent than those from any other 
gamma-ray observations. 



II. ANALYSIS PROCEDURE 

We begin our analysis by generating contour maps 
of the gamma-ray flux from the region surrounding the 
Galactic Center. These maps were generated using the 
latest (April 18, 2012) data release from the Fermi- 
LAT taken over the time period between August 4, 
2008 and June 19, 2012. Due to their superior point 
spread function, we use only front-converting events from 
the Pass 7 ultraclean class and, as recommended by 
the FGST collaboration, we include only events with 
zenith angles smaller than 100 degrees and do not in- 
clude events recorded while the Fermi satellite was tran- 
sitioning through the South Atlantic Anomaly or while 
the instrument was not in survey mode. 

Our raw maps are shown in the left frames of Fig. [I] for 
four different energy ranges between 300 MeV and 100 
GeV. In each map, ten contours are shown, distributed 
linearly between 2.45 x 10 -8 and 2.45 x 10 _7 cm _2 s _1 
sq deg" 1 (300-1000 MeV), 1.06 x 10~ 8 and 1.06 x 
10~ 7 cm- 2 s" 1 sq deg" 1 (1-3 GeV), 2.60 x 10~ 9 and 2.60x 
10- 8 cm- 2 s" 1 sq deg" 1 (3-10 GeV), and 3.60 x lO" 10 and 
3.60 x 10- 9 cm- 2 s- 1 sq deg- 1 (10-100 GeV). We have 
smoothed each of the maps at a scale of 0.5 degrees (the 
contour maps thus represent the average flux observed 
within a 0.5 degree radius of a given direction in the 
sky). 

In the right frames of Fig. [l] we show the maps af- 



ter subtracting the emission from known point sources 
and from the Galactic Disk, following the approach of 
Ref. [TTJ. In particular, we have subtracted a template 
map including all of the sources in the region contained 
within the Fermi Second Source Catalog [21] (adopt- 
ing central values for the intensity and location of each 
source, as reported in the catalog), with the exception 
of the bright central source, which cannot be easily dis- 
entangled from dark matter annihilation products in the 
case of a strongly cusped or contracted halo profile. To 
account for emission from the disk, we subtract a tem- 
plate with a morphology derived from the line-of-sight 
gas densities [351 US] as a function of galactic latitude 
(this template is only very mildly dependent on galactic 
longitude in the region of the sky being studied here). 
For details, see Ref. [TT]. As the morphologies of these 
subtracted backgrounds are not at all like those predicted 
from dark matter annihilation, we are in no danger of un- 
knowingly absorbing any would-be dark matter signal. 
Throughout this paper, dashed lines are used to denote 
negative contours of the same magnitude as assigned to 
the solid lines (resulting from oversubtraction) . 

After subtracting the contributions from the disk and 
from known point sources, central residuals remain in 
each energy range (outside of the inner ~2°, this sub- 
traction leaves only very modest residuals, typically on 
order of 10% or less of the residual flux in the innermost 
region). This central residual emission almost certainly 
includes some degree of contributions from the central su- 
permassive black hole [271 HB], unresolved point sources, 
and from cosmic ray interactions with gas (such as is 
observed at higher energies by HESS 29 , for example). 
To be conservative, we do not subtract any such com- 
ponents from this residual and derive all of our limits 
assuming that the entire residual flux could potentially 
be the products of annihilating dark matter. 

III. DARK MATTER ANNIHILATION IN THE 
INNER GALAXY 

The flux of gamma-rays in a direction, ip, from dark 
matter annihilations is given by: 

m = ^ d -^l Ami, a) 

where av and todm are the annihilation cross section and 
mass of the dark matter particle, respectively. dNj/dE 7 
is the gamma-ray spectrum produced per annihilation 
(as calculated using PYTHIA [30]), and the square of 
the dark matter density profile is integrated over the ob- 
served line-of-sight (los). This integral is often written in 
terms of the dimensionless function, J: 

J W-8^( o.3GeV/c m 3 )7o/ 2(0dL (2) 

In this section, we begin by considering a dark matter dis- 
tribution which follows the well-known NFW (Navarro- 
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FIG. 1: Contour maps of the gamma-ray flux from the region surrounding the Galactic Center, as observed by the Fermi 
Gamma-Ray Space Telescope. The left frames show the raw maps, while the right frames show the maps after subtracting the 
emission from the sources found in the Second Fermi Source Catalog (not including the central source) and from the Galactic 
Disk. All maps have been smoothed over a scale of 0.5 degrees. See text for more details. 



Frenk- White) profile 

p(r) oc 



where R s — 20 kpc is the scale radius of the halo. 
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And while the NFW profile is motivated by the results 
of numerical simulations [34], there is some (relatively 
mild) variation in the profiles that are predicted by such 
approaches. The results of the Via Lactea II simulation, 
for example, favor profiles with a somewhat steeper inner 



slope than NFW (p oc 



') |35j . while the Aquarius 



Project finds that the density slope varies with r [36] . 
The results of the Aquarius simulation have been used to 
motivate the commonly used Einasto profile: 



p(r) oc exp 



-2.0 

a 



1 



(4) 



where R s = 20 kpc is the scale radius of the halo. We 
will take a value of 0.17 for the parameter a. 

Neither the NFW nor the Einasto profile take into 
account the potentially important impact of baryonic 
physics on the distribution of dark matter in the inner 
Milky Way. Generally speaking, Milky Way-sized dark 
matter density profiles are expected to be contracted as 
a result of dissipating baryons, leading to the steepening 
of the inner profile [37]. The degree to which this effect 
is manifest depends on the fraction of the baryons that 
dissipate slowly by radiative cooling. 

Several state-of-the-art hydrodynamical simulations 
(performed by different groups, using different codes) 
have found that Milky Way-sized halos become signifi- 
cantly contracted by baryonic processes, increasing the 
density of dark matter in their inner volumes relative 
to that predicted by dark matter-only simulations (see 
Ref. [35] and references therein). These simulations, 
which include the effects of gas cooling, star formation, 
and stellar feedback, predict a degree of adiabatic con- 
traction which typically steepens the inner slopes of dark 
matter density profiles from 7 s» 1.0 to 7 ps 1.2 — 1.5 
within the inner ~10 kpc of Milky Way-like galax- 
ies [38j [39]. The resolution of such simulations is cur- 
rently limited to scales larger than ^100 parsecs [40]. In 
contrast, it has also been argued that strong feedback 
could result in the flattening of the inner slopes of galac- 
tic dark matter profiles [41j . This appears particularly 
likely in the case of low mass galaxies (ie. dwarf galax- 
ies), although it is less clear that such effects are likely 
to dominate in larger, Milky Way-sized, systems. 

In an effort to consider a wide range of possible bary- 
onic effects, we will derive limits using both contracted 
and cored halo profiles. To account for the possible ef- 
fects of baryonic contraction, we adopt a generalized form 
of the NFW profile: 



p(r) oc 



1 



(r/R s )-r{l + (r/R s )p-y)' 
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where 7 is the inner slope of the profile (7 = 1 recovers 
the NFW form) , and R s = 20 kpc is again the scale radius 
of the halo. In the following section, we will adopt values 
of 7 = 1.2 and 1.4, which we consider to represent mildly 
and significantly contracted profiles (although even larger 
values of the inner slope are sometimes found in hydrody- 



namical simulations [35], we will not consider such pro- 
files here). Alternatively, some simulations have found 
that resonant interactions between the stellar bar and 
the dark matter can flatten the density of the central 
cusp g2 i3] (see also, however, Refs. [riflg] ). With this 
in mind, we will also consider NFW-like profiles with flat, 
constant density cores within a radius of either 0.1 kpc 
or 1 kpc. 

Although profiles with inner slopes less steep than 
NFW (7 = 1) are sometimes found in hydronamical 
simulations with strong feedback, constant-density cores 
(7 = 0) are are rather extreme. Profiles with softened 
cusps (7 ~ 0.5 — 0.8, for example) exhibit much milder 
density suppression that is encapsulated in our constant- 
density cored profiles. Our results using these cored pro- 
files should thus be thought of as quite conservative, lead- 
ing to less bright dark matter annihilation signals than 
might otherwise be anticipated. 

At present, the distribution of dark matter in the 
inner Milky Way is not very strongly constrained by 
observations. In particular, the gravitational potential 
of the Milky Way's innermost kiloparsecs is dominated 
by baryons (stars and gas), which makes it difficult to 
strongly constrain the subdominant dark matter distri- 
bution in this region. In Fig. [2] we show the results 
of a recent study of microlensing and dynamical obser- 
vations of our galaxy [35] • While this study finds that 
these observational constraints are compatible with an 
NFW profile (7 = 1), they also allow for much steeper 
(7 ~ 1.8) or shallower (7 ~ 0.5) distributions. 

In each frame of Fig. [2] the best-fit profile is marked by 
a cross. To be conservative in the deriving of annihilation 
constraints, for each halo model we will normalize the 
dark matter density using the minimum (the lower 2a 
boundary) value found in Ref. 32J. These points are 
denoted by small diamonds in Fig. [2] 

In this study, we will derive dark matter annihilation 
constraints using six halo profile models. In the left frame 
of Fig. [3j we plot each of these density profiles (each nor- 
malized to the best-fit value as shown in Fig. [2]). Al- 
though these profiles are all very similar outside of the 
innermost few kiloparsecs of the Milky Way, they can 
vary considerably within the inner kiloparsec or so. In 
the right frame of this figure, we show the values of the 
density squared line-of-sight integral, J (see Eq. [2| , as a 
function of the angle away from the Galactic Center, tp, 
for each of these profiles. 

Upon first glance of this plot, one may expect the flux 
of gamma-rays from the innermost fraction of a degree 
around the Galactic Center to depend dramatically on 
the halo profile that is adopted. In reality, however, the 
finite angular resolution of Fermi (about 1° for a front- 
converting, 1 GeV photon) mitigates this variation to a 
significant extent. 
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FIG. 2: Constraints on the Milky Way's dark matter distribution from the combination of dynamical (rotation curves) and 
microlensing observations, as presented in Ref. 32 . The crosses in each frame denote the best-fit point. When deriving 
constraints on the dark matter annihilation cross section, we adopt halo profile normalizations which are at the lower 2a 
boundary of these constraints (shown as small diamonds). We thank the authors of Ref. [31] for providing the numerical values 
for the contours shown. 




FIG. 3: Left: Dark matter density profiles used in this study. Right: The value of the density squared line-of-sight integral, J, 
as a function of the angle from the Galactic Center, The contracted profiles are generalized NFW profiles with inner slopes 
of 1.2 and 1.4. The cored profiles are NFW profiles with constant densities within either 0.1 or 1.0 kpc. See text for more 
details. 
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IV. CONSTRAINTS ON THE DARK MATTER 
ANNIHILATION CROSS SECTION 

In this section, we derive constraints on the dark mat- 
ter's annihilation cross section for each of the halo profiles 
described in the previous section. To do this, we generate 
maps of the quantity J, convolved with the point spread 
function of Fermi, and again smoothed at the scale of 
0.5°. We then subtract these templates from the maps 
shown in the right frames of Fig. [I] with varying degrees 
of normalization. The results are shown in Figs. |4][9j In 
each case, the left and right frames include a contribu- 
tion from dark matter annihilations which is 40% less or 
more than in the center frame. 

As we increase the flux of dark matter annihilation 
products in a given energy range, eventually regions of 
the map begin to be oversubtracted (fluxes with negative 
values, as denoted by dashed contours). In the 300-1000 
MeV range shown in the top frames of Fig.|4j for example, 
we see significantly oversubtracted regions in both the 
middle and right frames. In either of these two maps, the 
oversubtracted regions are too large to be statistical fluc- 
tuations (no variations nearly this large are observed over 
the rest of the 40° x 40° region we have studied). From 
this procedure, we can determine the maximum flux of 
gamma-rays from dark matter annihilations in each en- 
ergy range, and convert this into an upper limit on the 
dark matter annihilation cross section, as a function of 
mass and dominant annihilation channel (quantitatively, 
we determine the upper limit as the point when we have 
oversubtracted a region that exceeds the scale of fluctu- 
ations along the disk at more than the 95% confidence 
level) . 

In Fig. [TUJ we show the resulting constraints under 
the assumption of either an NFW or Einasto profile. In 



Figs. 11 and 12 we show the constraints derived for con- 
tracted or cored distributions, respectively. In each case, 
we have conservatively adopted the minimum normaliza- 
tion for a given profile consistent with the Milky Way's 
rotation curve and microlensing constraints (see Fig. [5]). 
The constraints on the dark matter annihilation cross sec- 
tion would be about twice as stringent if we had adopted 
the central value for this normalization, as is often done. 

Although the strength of the resulting constraints de- 
pend on the dark matter distribution that is assumed, 
the variation from profile-to-profile is more modest than 
might be expected. For example, the constraints derived 
using a profile with a 1 kpc core are only a factor of a few 
weaker than those derived in the NFW case. The pres- 
ence of a 100 pc core has almost no effect on the resulting 
limits. The main reason for this is the limited angular 
resolution of Fermi (see the right frame of Fig.|3|. A sec- 
ond reason that limits the impact of the choice of halo 
profile is that the profiles which predict the highest an- 
nihilation rates (contracted profiles) also predict much of 
their signal to appear within the inner degree or so of the 
Galactic Center, where there is a significant gamma-ray 
flux observed. In contrast, profiles with a significant core 



Halo Model 


Energy Range 


f dE^^- (U.L.) 

m n M J 7 dE -< K ' 


NFW 
NFW 
NFW 
NFW 


0.3-1 GeV 
1-3 GeV 
3-10 GeV 
10-100 GeV 


l.OxlO -28 cm 3 s- 1 GeV _2 
3.6xl0 -29 cm 3 s" 1 GeV -2 
7.9xl0" 30 cm 3 s^GeV" 2 
9.5xl0" 31 cm 3 s" 1 GeV -2 


Einasto 
Einasto 
Einasto 
Einasto 


0.3-1 GeV 
1-3 GeV 
3-10 GeV 
10-100 GeV 


5.9xl0" 29 cm 3 s" 1 GeV" 2 
2.3xl0~ 29 cm 3 s" 1 GeV -2 
5.2xl0 -30 cm 3 s" 1 GeV -2 
5.2xl0~ 31 cm 3 s" 1 GeV -2 


Contracted (7=1.2) 
Contracted (7=1.2) 
Contracted (7=1.2) 
Contracted (7=1.2) 


0.3-1 GeV 
1-3 GeV 
3-10 GeV 
10-100 GeV 


4.6xl0 -29 cm 3 s" 1 GeV -2 
1.5xl0 -29 cm 3 s" 1 GeV -2 
2.8xl0" 30 cm 3 s^GeV -2 
3.4xl0~ 31 cm 3 s" 1 GeV -2 


Contracted (7=1.4) 
Contracted (7=1.4) 
Contracted (7=1.4) 
Contracted (7=1.4) 


0.3-1 GeV 
1-3 GeV 
3-10 GeV 
10-100 GeV 


9.0xl0" 30 cm 3 s^GeV" 2 
2.4xl0^ 30 cm 3 s^GeV" 2 
4.2xl0~ 31 cm 3 s" 1 GeV -2 
6.4xl0" 32 cm 3 s" 1 GeV -2 


Cored (J?c=100pc) 
Cored (7? c =100pc) 
Cored (7? c =100pc) 
Cored (7?c=100pc) 


0.3-1 GeV 
1-3 GeV 
3-10 GeV 
10-100 GeV 


l.lxlO" 28 cm 3 s" 1 GeV -2 
3.8xl0~ 29 cm 3 s" 1 GeV -2 
8.7xl0" 30 cm 3 s^GeV -2 
9.7xl0~ 31 cm 3 s" 1 GeV -2 


Cored (_R C =lkpc) 
Cored (i? c =lkpc) 
Cored (i? =lkpc) 
Cored (R c = lkpc) 


0.3-1 GeV 
1-3 GeV 
3-10 GeV 
10-100 GeV 


2.1xl0~ 28 cm 3 s" 1 GeV -2 
1.1 xlO -28 cm 3 s" 1 GeV -2 
3.0xl0 -29 cm 3 s- 1 GeV -2 
3.0xl0 -29 cm 3 s- 1 GeV _2 



TABLE I: The 95% confidence level upper limits on the quan- 
tity ((Tu/mn M ) J dE-ydN^/dE-y for various halo profiles and 
integrated over four different energy ranges. This table is 
intended to make it possible to derive limits for dark mat- 
ter models which are not explicitly considered in this pa- 
per, such as those in which the dark matter annihilates into 
combinations of different final states, or to any final states 
which are not considered here. For such models, one can 
use PYTHIA [30] to determine the spectrum of gamma-rays 
per annihilation, dN 1 /dE 1 , and then apply the constraints as 
presented in this table. In each case, we have conservatively 
normalized the halo profile to the minimum value capable of 
providing a good fit to the combination of the Milky Way's 
measured rotation curve and microlensing constraints |32| . 



are constrained primarily in the region 2-5° away from 
the Galactic Center, where the gamma-ray flux is lower. 

In order to present our results in a way that can be 
easily applied to dark matter models that are not ex- 
plicitly presented here (such as models which annihilate 
to combinations of different final states, or to final states 
that are not shown in our figures), we have present model 
independent limits in Table I. 

The constraints derived in this paper are somewhat 
more stringent than those presented in Ref. [IT]. This is 
largely due to the details of the dark matter distributions 
being considered. In Ref. [11] , halo profiles with a single 
power-law (p oc r -7 ) within the solar circle were adopted, 
recovering the NFW form in the R s — > oo limit. If R s is 



taken to have the conventional value used here (R s « 20 
kpc), the results of Ref. [TT] are in good agreement with 
those found in this study. 

Lastly, we note that some recent interest has been 
given to the possibility that the Milky Way's dark matter 
halo profile may peak at a location not coincident with 
the Galactic Center. This could be the case if our halo 
profile exhibits a large, nearly constant-density core |42j . 
In such a case, the resulting constraints would be only 
slightly different than those we have derived for the cored 
profile cases. On the other hand, it has been shown that 
the morphology of the 130 GeV line present within the 
FGST data d7J is best fit by a roughly Einasto-like dis- 
tribution, but centered around a point along the Galactic 
Plane approximately 1.5° away from the Galactic Center 
(/ = -1.5°) @H] (see also, however, Fig. 3 of Ref. 05]). 
In Fig. [l3j we consider this case and find a resulting con- 
straint of (cru/mQ M ) JiQQ e y dEydN^/dEy < 2 x 10~ 31 
cm 3 s -1 GeV -2 , which is more stringent than in the 
on-center Einasto case by a factor of approximately 2.5. 
This makes a dark matter annihilation explanations of 
the observed 130 GeV gamma-ray line very difficult to 
accommodate with dark matter candidates which pos- 
sess any significant annihilation cross section to non-line 
final states [5TM)2"] . 
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FIG. 4: Contour maps of the gamma-ray flux from the region surrounding the Galactic Center, after subtracting varying degrees 
of emission from dark matter distributed according to an NFW profile. As the flux of dark matter annihilation products is 
increased (moving from left-to-right), regions of the maps become increasingly oversubtracted (denoted by dashed contours). 
In this case of an NFW distribution, this occurs most noticeably in the regions approximately 1-2° north and south of the 
Galactic Center. 
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FIG. 7: Contour maps of the gamma-ray flux from the region surrounding the Galactic Center, after subtracting varying degrees 
of emission from dark matter with a significantly contracted (generalized NFW, 7 = 1.4) distribution. As the flux of dark 
matter annihilation products is increased (moving from left-to-right), regions of the maps become increasingly oversubtracted 
(denoted by dashed contours). 
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FIG. 8: Contour maps of the gamma-ray flux from the region surrounding the Galactic Center, after subtracting varying 
degrees of emission from dark matter distributed according to an NFW profile with a 100 parsec, constant-density core. As 
the flux of dark matter annihilation products is increased (moving from left-to-right), regions of the maps become increasingly 
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FIG. 9: Contour maps of the gamma-ray flux from the region surrounding the Galactic Center, after subtracting varying 
degrees of emission from dark matter distributed according to an NFW profile with a 1 kiloparsec, constant-density core. As 
the flux of dark matter annihilation products is increased (moving from left-to-right), regions of the maps become increasingly 
oversubtracted (denoted by dashed contours). 
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FIG. 10: The 95% confidence level upper limits on the dark matter annihilation cross section, for various annihilation channels, 
assuming a distribution which follows an NFW (upper frames) or Einasto (lower frames) halo profile. To be conservative, 
we have normalized the halo profile to the minimum value capable of providing a good fit to the combination of the Milky 
Way's measured rotation curve and microlensing constraints [32] (corresponding to a local density of p w 0.28 GeV/cm 3 or 0.25 
GeV/cm 3 in the upper and lower frames, respectively). For comparison, the horizontal line denotes the estimate for a simple 
thermal relic (av w 3 x 10" 26 cm 3 /s). 
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FIG. 11: The 95% confidence level upper limits on the dark matter annihilation cross section, for various annihilation channels, 
assuming a distribution which follows a mildly contracted (upper frames) or a significantly contracted (lower frames) halo 
profile. To be conservative, we have normalized the halo profile to the minimum value capable of providing a good fit to the 
combination of the Milky Way's measured rotation curve and microlensing constraints [31] (corresponding to a local density of 
p « 0.25 GeV/cm 3 or 0.22 GeV/cm 3 in the upper and lower frames, respectively). For comparison, the horizontal line denotes 
the estimate for a simple thermal relic (av w 3 x 10~ 26 cm 3 /s). 
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FIG. 12: The 95% confidence level upper limits on the dark matter annihilation cross section, for various annihilation channels, 
assuming a distribution which follows an NFW distribution with a constant-density 100 pc (upper frames) or 1 kpc (lower 
frames) radius core. To be conservative, we have normalized the halo profile to the minimum value capable of providing a 
good fit to the combination of the Milky Way's measured rotation curve and microlensing constraints |32] (corresponding to a 
local density of p ~ 0.28 GeV/cm 3 ). Notice that the limits derived for a profile with a 100 pc core are nearly indistinguishable 
from those derived in the NFW case. For comparison, the horizontal line denotes the estimate for a simple thermal relic 
(av « 3 x 10" 26 cm 3 /s). 
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FIG. 13: Contour maps of the gamma-ray flux from the region surrounding the Galactic Center, after subtracting varying 
degrees of emission from dark matter distributed according to a Einasto profile centered around the point (l,b) = (— 1.5°,0), 
as motivated by the morphology of Fermi's 130 GeV line [48| . As the flux of dark matter annihilation products is increased 
(moving from left-to-right), regions of the maps become increasingly oversubtracted (denoted by dashed contours). 
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V. COMPARISONS WITH OTHER DARK 
MATTER ANNIHILATION CONSTRAINTS 

In this section, we compare the constraints derived in 
this study to those from other gamma-ray observations. 
The most stringent of these constraints result from ob- 
servations of: 

• Dwarf Spheroidal Galaxies: Geringer-Sameth & 
Koushiappas pQ, and the Fermi Collaboration [5] 
have each published constraints on the dark mat- 
ter annihilation cross section from the stacked ob- 
servations of dwarf spheroidal galaxies. Both of 
these studies assume that the dark matter distri- 
bution in each dwarf is described by an NFW pro- 
file (in contrast, there is evidence from observa- 
tions [53H55] and simulations [56l [57] that favor 
cored profiles among dwarf galaxies) . Relaxing this 
assumption may weaken the resulting constraints, 
although only moderately . With this excep- 
tion, the constraints derived from dwarfs are robust 
and not subject to significant astrophysical uncer- 
tainties. In Fig. [14] we show these constraints, as 
derived from each of these groups [2 [2] (to present 
a fair comparison with our results, we have shown 
the constraint from Ref. [T] as derived using the 
minimum J factor, as we have normalized to the 
minimum acceptable halo density). 

• The Isotropic Gamma-Ray Background: The con- 
tribution to the isotropic gamma-ray background 
from cosmological dark matter annihilations could 
be significant. Depending on what one assumes 
about the quantity of substructure present in dark 
matter halos, constraints from the isotropic back- 
ground can be stringent. The uncertainties as- 
sociated with the issue of substructure are large, 
however. For example, in Fermi's analysis, results 
are presented for a range substructure models that 
produce annihilation constraints which vary by ap- 
proximately three orders of magnitude [S]. Only 
for the most optimistic substructure models are 
the constraints from the isotropic gamma-ray back- 
ground competitive with those presented in this 
study. In Fig. |14| we include the constraint from 
the isotropic background as presented in Ref. |8j as 
the "stringent constraint" for the case of substruc- 
ture model MSII-Sub 1. 

• Galaxy Clusters: Much like those derived from 
the isotropic gamma-ray background, dark matter 
annihilation constraints from galaxy clusters de- 
pend very strongly on the quantities of substruc- 
tures assumed to be present in such objects. The 
Fermi Collaboration's cluster analysis [5] was able 
to place only fairly weak constraints on the dark 
matter annihilation cross section (their most opti- 
mistic substructure model yielded av < 3 x 10 -25 
cm 3 /s for rriDM = 100 GeV and annihilating to bb), 



see Fig. [14] In contrast, the more recent analy- 
sis of Ref. [5] used a much more optimistic sub- 
structure model to produce constraints ~10 2 times 
more stringent. Given these enormous uncertain- 
ties, it can be somewhat difficult to interpret the 
constraints derived from clusters. 

In Fig. [14] we compare the constraints derived in this 
study to those derived from dwarf spheroidals, from the 
isotropic gamma-ray background, and from galaxy clus- 
ters. If we adopt an NFW halo profile (or an Einasto 
or contracted profile), the constraints derived from the 
Galactic Center are always the most stringent. Only if 
the dark matter halo profile of the Milky Way has a sig- 
nificant core (while dwarf galaxies retain their cusps) are 
constraints from dwarfs more stringent. The constraints 
from the Galactic Center are, for all dark matter masses, 
more stringent than those reliably extracted from the 
isotropic gamma-ray background or from galaxy clusters. 

VI. SUMMARY AND CONCLUSIONS 

In this article, we have used data from the Fermi 
Gamma-Ray Space Telescope to place constraints on the 
dark matter annihilation cross section. In doing so, we 
have considered a wide range of dark matter distribu- 
tions, each consistent with rotation curves and microlens- 
ing data. And while the strength of our constraints vary 
somewhat depending on which halo profile is adopted, 
we find that even in the most conservative cases (such 
as profiles with a significant, flat-density core) our con- 
straints are comparably stringent to those derived from 
observations of dwarf spheroidal galaxies (which were de- 
rived assuming NFW distributions). In all other cases, 
our constraints are significantly stronger than those de- 
rived from any other region of the sky. 

In an effort to be conservative in our analysis, we have 
chosen not to subtract a number of astrophysical compo- 
nents from the gamma-ray map of the inner Galaxy. In 
particular, we have not subtracted any component of the 
emission associated with the Galactic Ridge (as observed 
by HESS) or from the central gamma-ray point source (as 
observed by HESS, VERITAS, and MAGIC). Subtract- 
ing these components would have improved on the limits 
presented here by a factor of approximately two. Fur- 
thermore, we have in each case adopted normalizations 
for the dark matter halo profile which are the (2a) min- 
imum value consistent with dynamical and microlensing 
constraints. If we had instead used the central values for 
the overall dark matter density, our constraints would be 
more stringent by another factor of two. We have also not 
applied any boost factors or other enhancements due to 
dark matter substructure in deriving our results. The ex- 
istence of any such substructures would only strengthen 
our constraints further. 

Throughout this study, we have remained agnostic 
about the origin of the gamma-ray flux observed from 
the inner Galaxy. We only briefly mention here that 
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FIG. 14: A comparison of the upper limits on the dark matter annihilation cross section derived in this work to those from other 
gamma-ray observations. In particular, we show the constraints derived from the observations of dwarf spheroidal galaxies [HE], 
the isotropic gamma-ray background j8i, and from the Fornax galaxy cluster [3]. If we adopt an NFW halo profile (or an Einasto 
or contracted profile), the constraints derived from the Galactic Center are always the most stringent. Only if the dark matter 
halo profile of the Milky Way has a significant core (while dwarf galaxies retain their cusps) are constraints from dwarfs more 
stringent. The constraints from the Galactic Center are, for all dark matter masses, more stringent than those reliably extracted 
from the isotropic gamma-ray background or from galaxy clusters. 



the results presented here are in no way in conflict with 
those presented previously which find that annihilating 
dark matter can provide a good fit to the observed emis- 
sion pTTHT3| ITS] . In particular, Fermi's Galactic Cen- 
ter observations, coupled with observations of the Milky 
Way's radio filaments, are most easily explained by a 
dark matter particle with a mass of mdm ~ 7 — 10 GeV, 
an annihilation cross section of av ~ 5 x 10~ 27 cm 3 /s 
to charged leptons, and distributed in a somewhat con- 
tracted profile (p cx r -1 ' 3 ). 

Looking toward the future, we find very promising 
the possibility of the post-Fermi gamma-ray satellite, 
GAMMA-400 [59]. As GAMMA-400's overall effective 
area and acceptance will be comparable to that of Fermi, 
it will likely not be more sensitive to dark matter anni- 
hilations from flux-limited sources, such as dwarf galax- 
ies. With considerable improvements in both angular and 
energy resolution relative to Fermi, however, GAMMA- 
400 should be able to much better separate astrophysical 



backgrounds in the inner Galaxy from any dark mat- 
ter annihilation signal that is present. Furthermore, 
multi-wavelength studies of the Galactic Center, and 
progress from hydrodynamical simulations of dark mat- 
ter in Milky Way-like galaxies, could further strengthen 
the dark matter constraints that can be derived from the 
inner Galaxy. 
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